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We study about an evaporating process of black holes 
in S0(3) Einstein- Yang-Mills-Higgs system. We consider a 
massless scalar field which couple neither with the Yang-Mills 
field nor with the Higgs field surrounding the black hole. We 
discuss difi'erences in evaporating rate between a monopole 
black hole and a Reissner-Nortstrom (RN) black hole. Since 
a RN black hole is unstable below the point at which a 
monopole black hole emerges, it will transit into a monopole 
black hole as suggested via a catastrophe theory. We then 
conjecture the following: Starting from a Reissner-Nortstrom 
black hole, the mass decreases via the Hawking radiation and 
the black hole will reach a critical point. Then it transits to 
a monopole black hole. We find that the evaporation rate 
will increase continuously or discontinuously according to the 
type of phase transition, that is either second order or first 
order, respectively. After its transition, the evaporation will 
never stop because the Hawking temperature of a monopole 
black hole diverges at the zero horizon limit and overcomes 
the decrease of the transmission amplitude F. 

04.70.-S, 04.40.-b, 95.30.Tg. 97.60.Lf. 



I. INTRODUCTION 

For many years, there have been various efforts to find 
a theory of "everything" . One of the candidates is the 
superstring theory which has been the cause of much 
attention for last several years in the context of black 
hole thermodynamics. Since the discovery of black hole 
radiation by Hawking pi which is now called Hawking 
radiation, black hole thermodynamics takes its position 
beyond the analogy of ordinary thermodynamics. 

But Hawking radiation is a semiclassical phenomenon 
which means that space-time itself is treated classically 
and matter field is quantized around its metric. Although 
the gravitational field should also be quantized when the 
curvature radius gets as small as the Plancknian length 
{Ip '--^ 1.6 X 10~'^'^C7TT.), usually we ignore it and estimate 
the effect of Hawking radiation, e.g., as 7-ray sources of 
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the early universe 0]. When we consider the evapora- 
tion process of a Schwarzschild black hole, the Hawking 
temperature arises monotonically and Hawking radiation 
docs not stop, so classical physics will break down and 
quantum gravity effects should be considered. This is a 
serious unsolved problem which will be a key to quan- 
tum gravity. But if we think about Hawking radiation 
of a Reissner-Nortstrom (RN) black hole in the Einstein- 
Maxwell (EM) system, its fate is completely different be- 
cause its temperature will go down and the evaporation 
process will cease, if we assume that the electric charge 
is conserved. 

Such fates and related things have been investigated 
by many authors. But as for the black holes with non- 
Abelian hair ^^ , they have not been much investigated 
because their black hole solutions themselves are only 
obtained numerically, which takes much work compared 
with black holes analytically obtained. Another cause is 
perhaps due to their instability in which case the evapora- 
tion process need not be considered. But for some types 
of non-Abelian black holes, there exist stable solutions. 
Particularly, a monopole black hole which is found in 
SO (3) Einstein- Yang-Mill-Higgs (EYMH) system is im- 
portant in the context of Hawking radiation |[lO|-p^ . In 
the EYMH system, if we consider the evaporation pro- 
cess of the RN black hole, its fate will be rather different 
from that in the EM system, since it may experience a 
phase transition and become a monopole black hole. In 
contrary to the RN black hole, when a monopole black 
hole evaporates, the Hawking temperature rises mono- 
tonically like the Schwarzschild black hole and it may 
have the possibility to become a regular monopole. If 
this is the case, it may shed new light on the problem of 
the remnant of Hawking radiation. So we need to study 
its evaporation. 

This paper is organized as follows. In Sec. II, we intro- 
duce basic Ansatze and the field equations in the EYMH 
system. In Sec. HI, we briefly review black holes in 
the EYMH system and their thermodynamical proper- 
ties. In Sec. IV, we investigate the evaporating features 
of RN and monopole black holes in the EYMH system. 
In Sec. V, we make some concluding remarks and men- 
tion some discussion. Throughout this paper we use units 
c = h = 1. Notations and definitions such as Christof- 
fel symbols and curvature follow Misner-Thorne- Wheeler 



II. BASIC EQUATIONS 

We consider black hole spacetimc in the SO (3) EYMH 
model as 
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where k^ = SttG with G being Newton's gravitational 
constant. Lm is the Lagrangian density of matter fields 
which are written as 
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F^l^ is the field strength of SU(2) YM field and is ex- 
pressed by its potential AJJ as 

J^;. = d,Al - d.Al + ee'^^-AlAl , (3) 

with a gauge coupling constant e. <1>° is a real triplet 
Higgs field and Z?^ is the covariant derivative: 



A.*'' = 9„$'' + ee"''^^f,$' 



(4) 



The theoretical parameters v and A are a vacuum expec- 
tation value and a self-coupling constant of a Higgs field, 
respectively. To obtain black hole solutions, we assume 
that a space-time is static and spherically symmetric, in 
which case the metric is written as 
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where /(r) = 1 — 2Gm{r)/r. For the matter fields, we 
adopt the hedgehog ansatz given by 
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A^=^^e-VL^^, (a. = 1,2,3), (8) 

where r° and wf, are a unit radial vector in the internal 
space and a triad, respectively. 

Variation of the action (|l]) with the matter Lagrangian 
(PI) leads to the field equations 
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We have introduced the following dimensionless vari- 
ables: 



f = evr, rh — Gevm, 
and dimensionless parameters: 



vVG, \ = \/e' 
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Although the solution exists when v < Mpi, where Mpi 
is Planck mass, it can be described by a classical field 
configuration in the limit of a weak gauge coupling con- 
stant e, because its Compton wavelength ~ e/v is much 
smaller than the radius of the classical monopole solu- 
tion (~ 1/ev) in this case. Moreover, since the energy 
density is ~ e^w^ <C Mpi, we can treat this classically if 
we ignore the effect of gravity. The boundary conditions 
at spatial infinity are 

m(oo) ^ M < oo, 6{oo) — 0, h{oo) — 1, w{oo) — 0. 

(17) 

These conditions imply that space-time approaches a flat 
Minkowski space with a charged object. 

To obtain a black hole solution, we assume the exis- 
tence of a regular event horizon a,t r = rn- So the metric 
components are 



ruH = mirn) 



2G 



, Sh = S{rH) < oo. 



(18) 



We also require that no singularity exists outside the hori- 
zon, i.e., 
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For the matter fields, the square brackets in Eqs. (O) 
(O) must vanish at the horizon. Hence we find that 
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where 
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Hence, we should determine the values of wh, and hn 
iteratively so that the boundary conditions at infinity 
are fulfilled. 

Non-trivial solution does not necessarily exist for given 
physical parameters. However, for arbitrary values of v 
and A, there exists an RN black hole solution such as 
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M is the gravitational mass at spatial infinity and Q = 
2^v is the magnetic charge of the black hole. The ra- 
dius of the event horizon of the RN black hole is con- 
strained to be fii > Q. The equality implies an extreme 
solution. 

Around these black holes, we consider a neutral and 
massless scalar field which does not couple with the mat- 
ter fields, i.e., either Yang-Mills nor Higgs fields. This is 
described by the Klein-Gordon equation as 
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The energy emission rate of Hawking radiation is given 
by 
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where I and T[ijj) are the angular momentum and the 
transmission probability in a scattering problem for the 
scalar field <&. to and Th arc the energy of the particle 
and the Hawking temperature respectively. We define as 
S = -dM/dt. 

The Klein-Gordon equation (pif ) can be made separa- 
ble, and we should only solve the radial equation 



df*2 



X 



u? 



V' 



where 



V^ = 



f (l{l + l) 



fS' 



0, 



2(m — jfi'f) 



df 
df* 



fe 



(26) 

(27) 
(28) 



where ' denotes d/df and x is only the function of r. We 
need the normalization as Cj = oj/ev, V = V/ev. The 
transmission probability F can be calculated by solving 
the radial equation numerically under the boundary con- 
dition 
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where F is given as l/jAp. In our case, if we obtained the 
black hole solution, i.e., the shooting parameters wh and 
Hh, we should integrate ( p6| ) and (||)-(^ simultaneously. 



III. BLACK HOLES IN THE 
EINSTEIN- YANG-MILLS-HIGGS SYSTEM 

In this section, we briefly explain about black holes in 
the EYMH system, particularly about thermodynamical 
properties. We show the relation between horizon radius 
fn and gravitational mass M in Fig. 1(a). We denote 
RN and monopole black holes by a dotted line and a solid 
line, respectively. We chose as A = u = 0.1. In this figure, 
we can see that RN and monopole black holes emerge at 
the point B which does not change even if we change A. 
On the contrary, if we change iJ, the point B moves and 
disappears for large v (i,e., they do not emerge.). The 
precise are shown in |10|| . We only consider the parameter 
region where the point B exists, because it is shown in 
||ll| that the RN black hole becomes unstable via linear 
perturbation below the point B and later analysis showed 
that the RN black hole may transits to the monopole 
black hole only in this case. 

Before denoting the stability of monopole black hole, 
we show the magnification of Fig. 1 (a) around the point 
B in Fig. 1 (b). We find a cusp structure at the point 
A which exists for A < \crit- ^crit slowly depends on 
V. We can understand them via swallow tail catastro- 
phe Il3|. We can see in Fig. 1 (b) that for some mass 
range which corresponds to B to A, there appear three 
types of solutions (stable RN black hole, stable and un- 
stable monopole black holes) which suggests the violation 
of weak no hair conjecture. By contrast, for A > Xcrit, 
a cusp structure never appears and the monopole black 
hole solution merges with the RN black hole at the point 
B. From the analysis in [n3|, we can summerize the sta- 
bility of black holes as follows: 

(i) As for the RN black hole, the stability changes at 
the point B. It is stable or unstable according to be 
above or below the point B. 

(ii) As for the monopole black hole, when A < \crit, it 
is unstable along the curve from A to B, otherwise it is 
stable. When A > Xcrit, it is always stable. 

We also show the inverse temperature 1/Th in terms 
of the gravitational mass M/{M'^Jev) in Fig. 1(c). As- 
suming the conservation of charge and starting with the 
RN black hole, the point B is a key to the fate of the 
black hole. Because if we consider the RN black hole in 
the EM system, the RN black hole is always stable and 
the evaporation will cease at the extreme limit because 
the temperature vanishes there. While if we have the RN 
black hole in the EYMH system, the RN black hole be- 
comes unstable below the critical point B and will change 
to a monopole black hole by either second- or first- order 
transition according to A > Xcrit or A < Xcrit- After this 
transition, because the temperature diverges to infinity 
at the rn ^ ^ limit, so we may not stop evaporating a 
monopole black hole and find that one of the candidates 
for the remnant is a self-gravitating monopole. We show 
the diagram for A = u = 0.1. But the results are qual- 
itatively same for other parameters (see Fig. 9 (a) in 



The criterion stable or unstable is based on casastrophe 
which coincide with the analysis via linear perturbation. 
So if we think evaporation process and time evolution of 
them, the results may change. But it may be laborious to 
calculate such evolution, here we consider mainly about 
transmission amplitude of a scalar field assuming black 
holes as a background. 



IV. EVAPORATION OF BLACK HOLES IN THE 
EINSTEIN- YANG-MILLS-HIGGS SYSTEM 

Even if we assume that the background spacetime is 
fixed and ignore backreaction, it is not evident to pre- 
dict the final fate of black holes. Naively speaking, the 
temperature is the main cause to decide evaporation pro- 
cess. But the transmission amplitude T may also affect 
the results. For example, a dilatonic black hole in the 
EM-dilaton system has different properties via a coupling 
constant a of the dilaton field to the matter field [Esj . If 
a > 1, the temperature of the black hole diverges and 
the effective potential V'^ grows infinitely high simulta- 
neously at the extreme limit [^. In this case, it is not 
evident how to decide whether or not the emission rate 
diverges. In |I^, it turned out that the divergence of 
the temperature at the extreme limit overcomes that of 
the effective potential, resulting in a divergence of the 
emission rate. 

In the case of a monopole black hole, it is not even evi- 
dent whether or not the effective potential diverges at the 
rjy — > limit, because its solution is only obtained nu- 
merically. Fig. 2 shows the effective potential (V/ev)'^ in 
terms of the radial coordinate for v/Mpi = 0.05, A/e^ = 1 
and rn /{ev)~^ = 0.1, 0.3. Taking the horizon radius as 
smaller, the potential V'^ becomes larger and our numeri- 
cal results suggest that its potential diverges within that 
limit, so we must analyze the emission rate to decide 
whether or not the evaporation will stop. Another in- 
teresting point is how Hawking radiation changes at the 
point B in the transition process of an RN black hole to 
a monopole black hole. Near the point B we can suggest 
something concrete under some assumptions. For this, 
we assume the following: (i) The coupling constant e is 
small enough so that we can treat the gravitational field 
classically at the point B. (ii) A discharge process does 
not occur during the evaporation fls] . (iii) The coupling 
of the matter fields (YM and Higgs fields) to the scalar 
field even if it exists, the results would not be much af- 
fected. The last ansatz may seem to be strong, but it 
might turn out to be true near the point B because the 
monopole black hole around there is very close to the RN 
black hole, which is a vacuum solution. In fact the field 
strength of YM field for such a monopole black hole is 
much smaller than that in the other cases. 

Before seeing such properties, we show the trans- 
mission probability F of a RN black hole and a 



monopole black hole in terms of lo for the horizon ra- 
dius TH/iev)-^ = 0.55 and v/Mj,i = 0.05, A/e^ = 0.1, 
1 in Fig. 3 (a). We depict only I — 0, 1 modes. The 
RN black hole, monopole black hole with A/e^ = 0.1, 
monopole black hole with A/e^ — 1 correspond to the 
dotted lines, the solid lines and the dot-dashed lines, re- 
spectively. As we see the dominant contribution for the 
Hawking radiation is Z = 0, because the contributions 
of the higher modes are suppressed by the centrifugal 
barrier. In what follows, we will ignore the contribu- 
tions from I > 2. Although F is the largest for an RN 
black hole, these are almost indistinguishable. So, when 
evaluating the value of the emission rate of a monopole 
black hole, we may conclude that F may not be the main 
origin of the difference from that of the RN black hole. 
However, for a monopole black hole with a smaller hori- 
zon radius, its difference from an RN black hole becomes 
clear. Fig. 3 (b) shows the same diagram in Fig. 3 (a) 
with the same parameters v and A but with a smaller 
horizon radius rn ~ 0.3/ev. We can see the difference 
clearly. It is because the size of non-trivial structure be- 
comes larger compared with the horizon radius for the 
monopole black hole of smaller horizon. We also show 
the examples of the radiation spectrum —cPM/dtdio in 
terms of w in Fig. 4 (a), (b). The parameters correspond 
to Fig. 3 (a), (b). One can see the reason why we can 
ignore the contribution for Z > 2. Actually, the difference 
caused by it is below 1% from one of our calculation. 

We return to the first concern, i.e., what happens when 
the horizon radius changes via Hawking radiation. In Fig. 
5 (a), we show the emission rate S in terms of the grav- 
itational mass M of an RN black hole and a monopole 
black hole for v/Mpi = 0.05, A/e^ == 1. The difference 
between an RN black hole and a monopole black hole 
may be caused by the Hawking temperature Ty, because 
the emission rate oc Tfj. In Fig. 5 (b), we show the 
emission rate S in terms of the horizon radius fn for the 
same parameters in Fig. 5 (a). This diagram strongly 
suggest that the evaporation will not stop even at the 
r/f — > limit. This resembles the situation of the dila- 
tonic black hole at the extremal limit, i.e., whether or not 
evaporation stop depends only on Th |1^. Thus when- 
ever we think black hole space-times as background, the 
transmission amplitude does not change the scenario es- 
timating from the temperature. 

We estimate the time scale of the evaporation using 
Ttime = M/S as the indicator, and show this time scale 
in terms of the gravitational mass M in Fig. 6 for the 
same solutions in Fig. 5 in the CGS units. Near the 
bifurcation point B, Tume ^ lO^'^^'/e'^ second^ and below 
the bifurcation point, the time evolution of the monopole 
black hole is completely different from that of an RN 
black hole. 

Next, we study the evaporation rate near the bifurca- 
tion point when either the first or second order transition 
to the monopole black hole occurs. In Fig. 7, we show 
the emission rate S in terms of the gravitational mass 
M of RN and monopole black holes for v = 0.05Mp/, 



A/e^ = 0.1, 0.2, 0.3, 0.4, 1 near the bifurcation point. 
The curves from B to A correspond to the emission rate 
of the unstable branch. In the case that the transition is 
first order, the emission rate will also change discontinu- 
ously as shown by an arrow in Fig. 7. 

It may be interesting to ask in which direction the RN 
black hole jumps to the monopole black hole or how much 
black hole entropy (or radius of event horizon) will in- 
crease after the phase transition. In order to analyze this 
problem properly, we have to include a back reaction ef- 
fect of Hawking radiation, which is very difficult and has 
not yet been solved. However, we may find some con- 
straints through the following considerations. Because 
the RN black hole emits particles, it will lose some of 
the gravitational mass. But, whether the horizon radius 
increases or not may depend on two time scales, i.e. the 
evaporation time and the transition time. If we apply 
the catastrophe theory, the entropy of the black hole will 
increase and the horizon radius will increase. But this 
analysis is based on the assumption that the change of 
black hole states can be treated adiabatically. Since the 
coupling constant e is so small that each state can be 
described by a quasi-stationary solution, we may expect 
that the horizon radius will increase after the transition. 

We also confirm that an other choice of values of v and 
A does not provide any serious difference in the evapora- 
tion process. In Fig. 8, we show one of the interesting 
cases (A/e^ = 0.1 and v = 0.2Mpi) where the bifurcation 
point B is very near the extreme RN solution. In this 
case, this diagram suggests that the RN black hole first 
almost ceases the evaporation process and becomes close 
to the extreme one, and then it transits to a monopole 
black hole and will start to evaporate again. In other 
parameters, these diagrams are basically similar to the 
above cases: 

(i) The emission rate diverges at the fn—'O limit. 

(ii) Near the bifurcation point B, the emission rate 
changes continuously or discontinuously according to A 
is above or below Xcrit- 



V. CONCLUSION AND DISCUSSION 



that the evaporation rate changes continuously or dis- 
continuously depending on whether the transition which 
occurs is second- or first-order. Our results suggest that 
the Hawking radiation near the bifurcation point is de- 
termined only by the temperature of the black hole. We 
can understand this as follows. Because in particular, 
around this region, we find little difference in the trans- 
mission probability between a monopole black hole and 
an RN black hole. 

(ii) When the horizon radius becomes small, the trans- 
mission probability of a monopole black hole becomes 
small compared with that of an RN black hole. Though 
it can not stop evaporating in our analysis because the 
increase of temperature of a monopole black hole at the 
rn —^ limit, quantum effects of gravity may cause a 
serious effect on it and would overcomes the decrease of 

r. 

We finally remark on some subjects which we leave to 
the future. When we consider the fate of an RN black 
hole via Hawking radiation, we may take into account the 
effects of charge loss if it is to be expected, and have to 
include a coupling to the YM field or Higgs field before we 
consider the effects of quantum gravity. The second is the 
concern with the critical behavior ||l^ . There are several 
works about it in the EYM or Einstein-Skyrme systems 
P0| , pi| , in which the Schwarzschild black hole is the most 
stable one. But in the EYMH system, since a monopole 
black hole becomes more stable than the RN black hole 
below a certain critical mass, it would be interesting to 
study the critical behavior in the EYMH system. Finally, 
it may be more interesting to look for the "real" critical 
behaviour in our present phase transition via Hawking 
evaporation. Those are under investigation. 
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We have considered an evaporation process of the RN 
and monopole black holes in the EYMH system. We 
have analyzed a real massless scalar field which couple to 
neither the NA field nor the Higgs field. We may suggest 
how RN and monopole black holes evolve through an 
evaporation process in the EYMH system. We have the 
following results. 

(i) We investigated the evaporation process of an RN 
black hole, in particular, near the bifurcation point where 
this merges with a branch of monopole black holes. Since 
the RN black hole becomes unstable there, we expect 
that it transits into a monopole black hole. This transi- 
tion will be first- or second-order according to whether A 
is smaller than some critical value Xcrit or not. We show 
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FIG. 1. (a) the gravitational mass M/{Mpi/ev) and the 
horizon radius rH/iev)^^ relation and (b) its magnification 
diagram near the point B and (b) the inverse Hawking tem- 
perature {Tn/eAIpi)'^ in terms of the gravitational mass 
M/{Mpi/ev) of the monopole black hole with A/e^ = 0.1 
(solid lines) and of the RN black hole (dotted lines). We 
choose v/Mpi = 0.1 in these diagrams. At the point B, the RN 
black hole becomes unstable and change to monopole black 
hole. This process is first or second order correspond to A 
is larger than Xcrit or not. Fig. (c) shows that temperature 
of the monopole black hole diverges at the rn ^> limit like 
the Schwarzschild one. Though we showed the one example, 
which is irrelevant to the results (see Fig. 9 (a) in |13].). 
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FIG. 2. We show the potential 1^^/(611)^ in terms of 



the radial coordinate with 

-1 



\/e' = 1, 



v/Mpi = 0.05 and 
r_tf/(e«)~^ — 0.1, 0.3 for I — mode. It shows that when 
horizon radius th becomes small, potential barrier becomes 
large. It seems that if we take th — > 0, potential barrier 
diverges. 





FIG. 3. Energy dependence of the transmission amplitude 
r for monopole black holes with A/e^ = 0.1 (solid lines) and 
A/e^ = 1 (dot-dashed lines) and of RN black holes (dotted 
lines) for / = 0, 1 mode which make the main contributions 
to the Hawking radiation. We choose v/Mpi — 0.05 in these 
diagrams and (a) rn/iev)^^ — 0.55 (b) r_ff/(e«)~^ — 0.3. 
Though the RN black hole has the largest F among them, 
these are almost indistinguishable in (a). But when horizon 
radius becomes small, their difference becomes large because 
YM field and Higgs field outside the horizon becomes large 
and contribute to the black hole structure. 
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FIG. 4. Energy spectrum of the —d m/dtduj for black holes 
corresponding to Fig. 3 (a), (b). They shows that why we 
can neglect the contribution from I > 2. 
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FIG. 5. The emission rate H/ew (a) in terms of the gravi- 
tational mass M/{Mpi/ev) (b) in terms of the horizon radius 
rn/iev)^^ for RN and monopole black holes for v/Mpi — 0.05, 
A/e^ = 1. If we assume the charge is conserved, though the 
RN black hole will stop evaporating at the extreme limit, 
monopole black hole does not stop evaporating as in the 
Schwarzschild black hole. If we consider the effects of quan- 
tum gravity, the results near rn ~ Ip may be changed. 
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FIG. 6. The evaporation time scale Tume in terms of the 
gravitational mass M in the CGS units for the same solutions 
in Fig. 5. 
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FIG. 7. The emission rate H/e« in terms of the gravita- 



for RN and monopole black holes 
= 0.1, 0.2, 0.3, 0.4, 1 near the bi- 



tional mass M/^MlJe 

for v/Mpi = 0.05, A/e 

furcation point B. The lines from _B to A correspond to the 

emission rate of the monopole black hole which is thought to 

be unstable. So when the transition from the RN black hole 

to the monopole black hole occurs, emission rate will raise 

above the points A. 




FIG. 8. The emission rate H/e« in terms of the grav- 
itational mass M/{Mli/ev) for v/Mpi = 0.2. We chose 
A/e^ = 0.1. In this figure, we considered the situation bifurca- 
tion point almost coincide with extreme RN black hole. This 
figure shows that evaporation feature drastically changes near 
the extreme point because of the transition into the monopole 
black hole. 



